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Abstract 

One of the unsolved problems in cosmic ray (CR) physics is the small radial gradient of the 7-ray intensity 
compared to the inferred CR source distribution in the Galactic disk. In diffusive CR propagation models the 
most natural explanation is very efficient spatial mixing due to MHD turbulence in the interstellar medium. 
However, even in the most favorable case of a very large diffusive CR halo the 7-ray gradient is still larger 
than observed. In our view the small 7-ray gradient could be the result of strong advection by a Galactic wind. 
We show that a small 7-ray gradient can be obtained, if the diffusion region does not extend far above the 
Galactic plane. Important ingredients of our model are: (i) anisotropic CR diffusion, (ii) strong radial and 
<Q vertical gradients of the advection velocity (due to faster winds above higher CR source density regions). 

1 Introduction: 

The distribution of the 7-ray emissivity in the galactic disk bears important information on the CR origin, 
because Galactic diffuse 7-rays result from interactions of CRs with interstellar gas. Thus, this distribution 
depends critically on the conditions of CR propagation in the Galaxy. For 7-ray energies above 100 MeV, 
the main production process is 7r°-decay, resulting from nuclear collisions between high energy particles and 
interstellar matter. The study of diffuse 7-ray emission has shown that the nucleonic component of cosmic 
rays is more or less homogeneously distributed over the entire Galactic disk. If the spatial distribution of 
CRs were uniform, the 7-ray emission should map perfectly well the distribution of interstellar hydrogen. 
However, the distribution of the cosmic ray (CR) sources (most likely supernova remnants (SNRs) for particle 
energies below 10 15 eV) is far from homogeneous. Taken at face-value, the discrepancy must arise during 
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2 ■ the propagation of CRs from their sources through the interstellar medium. The difficulty in interpreting the 
data correctly is due to the fact, that the coefficients of CR transport are not well known, and therefore models 
(plus their respective assumptions) are used to bridge the gap between the measured 7-ray parameters and the 
unknown parameters of CR origin. 

The most common models for the interpretation of the spatial CR distribution in the disk are phenomeno- 
logical in nature and based on CR diffusion. It is assumed that particle propagation in the interstellar medium 
can be described by an isotropic diffusion process due to CR scattering off magnetic field fluctuations. The 
value of the diffusion coefficient (or tensor) is estimateded from observational data. In many cases this model 
gives an acceptable interpretation of experimental results (see e.g. Berezinsky et al., 1990). 

The interpretation of the 7-ray gradient in the framework of an isotropic diffusion model, in which the 
CRs are injected by a radially non-uniform SNR distribution (e.g. Case & Bhattacharya, 1996) leads to the 
following conclusion: in order to obtain a rather uniform CR distribution in the disk one has to assume that 
the Galaxy is surrounded by a huge halo whose vertical extension is larger than 10 kpc (Bloemen et al., 1993). 
Only then diffusion mixes CRs efficiently enough, so that their distribution in the disk differs appreciably from 
that of their sources. The numerical calculations (see Dogiel & Uryson, 1988) show however, that even in the 
case of a very extended diffusion halo the calculated gradient of CRs is larger than derived from COS-B and 
EGRET data. Therefore, if we believe that the observations of the SNR and 7-ray distributions are correct, we 
should seek alternative explanations for the gradient data. 



2 Effects of CR Transport by a Galactic Wind: 



A completely different scenario of CR escape into intergalactic space is implied by the so-called galactic 
wind models. It is shown that the combined pressure of thermal plasma, CRs, magnetic fields and MHD waves 
can lead to a secular escape of gas and CRs in normal spiral galaxies, provided that the coupling between 
scattering waves and energetic nucleons is strong. In this case there is a net forward momentum transfer from 
CRs to the gas via the frozen-in waves as a mediator (see Breitschwerdt et al., 1991; 1993). In wind models 
the coupling between MHD-waves and CRs is generally assumed to be strong everywhere in space except for 
a rather narrow region surrounding the galactic plane where CRs propagate predominantly by diffusion. In 
the wind CRs are transported along with the gas and the waves to large distances from the disk. As numerical 
simulations have shown, these models exhibit two effects that may provide the long-sought explanation for the 
7-ray gradient variation problem, both of which are not present in static halo models. The first one is that the 
diffusion tensor may be anisotropic. The reason is that the component along the magnetic field lines depends 
on the CR energy, as Dm oc E om , and the component perpendicular to the lines D± is energy independent. As 
a result in the case of a strong galactic wind the location, z c , of the halo boundary depends on the energy as: 
z c oc E 0M and the local CR life time T cr ~ ^ oc E 06 . The diffusion component Dm rises with energy so the 
diffusion is one dimensional at high energies. 

The second, and more important, new effect is the radial variation of the wind velocity (or mass flux, if the 
gas density is constant) of the galactic wind along the galactic plane, which may give rise to an almost uniform 
CR distribution in the disk, even for a strongly non-uniform CR source distribution, provided that the diffusion 
halo height is small. The reason is that, depending on the local CR source strength (i.e. the number of SNRs 
per unit volume and hence the number density of CRs at a given energy, n(r, z, E)), the CR energy density, e c 
or its pressure Pq, will vary accordingly with galactocentric radius. Thus an increase in e c will also increase 
the galactic wind mass flux or wind velocity, respectively, and hence decrease the distance to the advection 
boundary, given by D(R C , z c , E)/ (V(R C , z c )z c ) ~ 1, where R and V = u + Va denote the halo radius and 
the CR transport speed with respect to an Eulerian frame of reference. Consequently, the CR storage volume 
will be reduced, thus facilitating CR escape locally and thereby decreasing the number of nuclear collisions, 
which produce 7-rays via 7r°-decay. 

3 Uniform Disk Distribution of CRs in the Wind Model: 

To illustrate the behaviour of the solution of the advection-diffusion equation we discuss two simple cases. 
The first one describes one-dimensional diffusion with the advection velocity V depending on the radial coor- 
dinate r only, as V(r) = Vq ■ f\ (r), with Vq being a constant and f\ (r) an arbitrary function of r. The equation 
to solve for the CR distribution function N(r, z) is 



with the boundary condition N = at z c = 1. Here Q is the non-radial part of the CR source distribution, and 
/2(r) is also an arbitrary function of r, which describes the radial SNR distribution in the disk. The solution 
of Eq. Q is 



We see that for the case of weak advection Vq <C 1 the solution is N (5/2(1 — z)/D, and the CR gradient 
depends on the source distribution / 2 only. However if Vq S> 1 then 



and in this case the gradient is a strong function of the velocity distribution too. In the special case /i(r) = 
/2(r) we obtain N = const. 
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Let us investigate a more complex case in which the velocity is also a function of z. For a velocity distri- 
bution V(z,r) = 3 Vq z fi(r), with some suitable constant Vq, the solution of the one-dimensional diffusion 
equation was obtained by Bloemen et al., 1993. For strong advection the solution converges to 

N(z = 0, r) ~ A Qh{r \ - ; (4) 

we retrieve N = const, if /i(r) = /f ( r )> where A is a constant. 

The diffusion coefficient D is assumed to be constant, although we cannot exclude its spatial variation, 
which can change the CR density distribution in the disk as well. 

4 Analytical Solution of the Diffusion- Advection Equation: 

We now want to work out a general solution of the two-dimensional CR transport equation for nucleons, 
which reads 



V {fi(x, E)VN - u(x)N) - A Q W(x) En]--^(^n)=Q(E,x 



(5) 



where x denotes the spatial coordinates, and the diffusion tensor D_ in cylindrical coordinates with axial sym- 
metry is given by 

£ = {l>Z dZ) = { D,°E°)- (6) 

Our prime interests here are the anisotropic diffusion and advection, rather than the energy dependence and 
so, for convenience, we set a = 0. For the nucleon component other than adiabatic losses are negligible 
(dE/dt = 0); thus Eq. $5^ in axial symmetry becomes 
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where we have chosen a spatial variation for the velocity field as 

**> =(£)=( ) • (8) 

and we have used a power law spectrum for particle injection by the disk sources. The choice of u{x) seems 
rather arbitrary at first glance, and it can indeed be quite a complicated function of r and z. However, numerical 
calculations of galactic winds (Breitschwerdt et al., 1991) show a similar spatial variation of the velocity field 
like the one we have chosen. 

The full analytical solution od Eq. © is rather cumbersome, and will be discussed in detail elsewhere 
(Breitschwerdt et al., 1999). Here we present the result for the CR distribution function in the disk (z = 0), 
subject to the boundary conditions N(r, z = ±oo) = and N(r = ±oo, z) = 0: 
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Here, (3 n = -2n - 3A - 2- s /A(j + 6re), % = ±(1 - 4n - 3A - 2 V / A{~/ + 6n)), 7™ = |(2(2m + 1) - 
3,4 + 1 + 2^^(7-3(2771 + 1))), K n = -2n - ^(7 + 6n), dK n / da = 1 + 3/V(2y /A(7 + 6nJ), 
a m = 2m + 1 - 3A + 2^4( 7 - 3(2m + 1)), iv~ m = 2m + 1 + v^(7 - 3(2m + 1)), dK m /d/3 = 
1 - 3/1/(2t/A(7 - 3(2m + 1))), where n,m= 0,1,2,..., D z = D±, A = V /k t , B = \J D z / K r , f = Br, 
and 7 = 70 + 2. 

It is instructive to look at the asymptotic behaviour of the solution at different radii. If the CR sources 
occupy a limited volume of the disk bounded by a radius a, and we let Q(r') = Qr'~ q <d(a — r') (where q 
is a fit to the observed SNR distribution, and denotes the step function), then for r > a the function N is 
completely determined by the second integral which is a constant. Indeed from contour integration we find 
that the first integral / r + °° = and the second integral J r = Jq = const for r > a. Then from a simple 

analysis we see that N(r) oc r -1- "^/^ 7-3 ) and N(r) oc r _1 if A 1. This is just what is expected for purely 
diffusive particle propagation. In the case of strong advection, i.e. A S> 1, we have N(r) oc r^V^ -3 ); such 
a strong drop of CR density far away from the sources is due to adiabatic energy losses, since, if we neglect 
them, 7 = 3 formally, and N(r) oc r _1 , almost independent of the value of A. 

In the vicinity of the source region (r <C a) the function N is determined by the two integrals (cf. Eq. Q) 
which can be written as 

N(r) oc r ^-^) dr > _ r -i-^(7-3) f/^+^^V . (10) 
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For the case of strong advection i» 1 the first integral is determined by the lower limit and the second one 
by the upper limit, independent of the source distribution. Then we have N(r) oc r 2 ~ q and the CR density is 
almost in the disk region close to the sources if q = 2. 

If A <C 1 then the values of the integrals are determined by the source distribution. If the CR sources 
are concentrated towards the centre of the disk, q > 2, then the CR distribution is determined by the second 

integral and N(r) oc r -1- "^^ 7-3 ). If the sources are uniformly distributed in the disk (q = 0) then N oc 
(a 2 — r 2 ). Thus from these analyses we see that a more or less uniform CR distribution can be expected in the 
disk, if advection is strong, i.e. A S> 1. 

5 Conclusion: 

If the advection velocity is proportional to the CR source density in the disk, then the propagation char- 
acteristics of the CRs at each point of the disk are strongly determined by local conditions. In this case we 
can obtain a radially uniform CR distribution, if the halo boundary is close to the disk. Therefore local CR 
characteristics should vary strongly from point to point in the disk. Moreover, any attempt to estimate global 
parameters of CRs (such as the halo height or diffusion coefficient) from distribution of their density in the 
disk is misleading, since mixing of CRs in the Galaxy may be rather weak. Therefore one may conclude that 
in the case of strong CR advection, the halo extension derived from CR nuclear data reflects only a local halo 
extension near Earth and the value derived from 7-ray data and pure diffusion models may be an artifact. 
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